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ABSTRACT 

We report the dramatic mid-infrared brightening between 2004 and 2006 of HOPS 383, a deeply 
embedded protostar adjacent to NGC 1977 in Orion. By 2008, the source became a factor of 35 
brighter at 24 /im with a brightness increase also apparent at 4.5 /im. The outburst is also detected 
in the submillimeter by comparing APEX/SABOCA to SCUBA data, and a scattered-light nebula 
appeared in NEWFIRM K s imaging. The post-outburst spectral energy distribution indicates a Class 
0 source with a dense envelope and a luminosity between 6 and 14 Lq. Post-outburst time-series mid- 
and far-infrared photometry shows no long-term fading and variability at the 18% level between 2009 
and 2012. HOPS 383 is the first outbursting Class 0 object discovered, pointing to the importance of 
episodic accretion at early stages in the star formation process. Its dramatic rise and lack of fading 
over a six-year period hint that it may be similar to FU Ori outbursts, although the luminosity appears 
to be significantly smaller than the canonical luminosities of such objects. 

Subject headings: Stars: formation - - Stars: protostars — circumstellar matter — Infrared: stars 


1. INTRODUCTION 

Early infrared studies of large samples of young stars 
revealed that p rotostars were under-luminous compared 
to predictions l)Kenvon et al.l ll990l) . In order to form 
a solar-mass star in 10° yr, the time-averaged mass ac¬ 
cretion rate onto the star over this period needs to be 
10 -5 M q yr, implying total luminosities in excess of 
10 Lq. Median protostellar lumino s ities, however, are 
of order 1 L ,n dKenvon et al.l 119901 lEvans et al.1 120091 
iKrvukova et al.ll2012f) . One means of resolving the dis¬ 
crepancy is episodic accretion , in which the luminosity of 
a forming star is usually ~ 1 Lq as observed, but a series 
of relatively brief, dramatic spikes in the accretion rate 
and luminosity over the star-formation peri od su pply the 
req uisite mass for a sun-like star (lHartmann &; Kenvonl 
fl996h . 

During the 20th century, several young stars were ob¬ 
served to undergo bursts consistent with the episodic 
accreti on hypothesis, beginning with FU Ori in 1936 
(jWachmanni [ 1939 1. With extensive optical to infrared 
surveys of star-forming regions by, e.g., the Palomar 
Transient Factory and the Spitzer Space Telescope, and 
with careful work by amateur astronomers such as R. 
Persson (V733 Cep), J. McNeil (V1647 Ori), and A. 
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Jones (EX Lup), many additional outbursts have been 
discovered ov e r the past decade. The recent review of 
lAudard et al.l (l2014f) catalogs 26 eruptive young stars. 

While the first outbursts to be observed resembled 
Class II young stellar objects (YSOs), with luminous 
accretion disks but weak to absent circumstellar en¬ 
velopes, further discoveries have extended the outburst 
phenomenon to the envelope-embedded Class I proto¬ 
stars. Such protostella r outburst s include V346 Nor 


dGraham k. Krone] 1 1 9851) . OO Ser dHodapp et all 


iKosnal et al.ll2007F V1647 Ori (jReinurth fc As nil'll 


1996 


2004: 


lAbraham et al.112004 iBri ceno et akl I2004D. V277 5 Ori 


(|Caratti o Garatti et al.l 12011b Fischer et al.l 120121) . and 


V900 Mon I Reimirth et al.l 120121) . By obtaining 70 /im 
light curves for protostars in the Orion Nebula Cluster, 
iBillot et~akl (|2012T ) also found evidence for variable ac¬ 
cretion below the level associated with bursts, but still 
in excess of the 10% level, in eight of 17 protostars. 

The relative importance of episodic ac cretion to the 
star-formation process is sti ll under debate (|Baraffe et al.l 
120091 [Hosokawa ct al. 2011). Is the majority of the mass 
of a typical solar-mass star accreted in stochastic out¬ 
bursts or in a smooth process of infall from a circum¬ 
stellar envelope, through an accretion disk, and onto the 
star? The likelihood of the former scenario increases as 
new outbursts are discovered at earlier stages of the star 
formation process. Here we announce the discovery of 
an outburst in HOPS 383, a Class 0 protostar with the 
reddest spectral energy distribution (SED) yet observed 
for an outburst, implying that it may be the youngest 
known episodic accretor. 


2. OBSERVATIONS 


In Orion, we performed a search for variability consis¬ 
tent with episodic accretion by comparing photometry 
of protostars at 3.6, 4. 5, and 24 pm f rom t he 2004- 
2005 Spitzer survey by iMeeeath et al.l (120121 ) to pho¬ 
tometry at 3.4, 4.6, and 22 /im from the 2010 sur¬ 
vey of the Wide-Field Infrared Survey Explorer ( WISE ; 
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IWright et al.ll2010h . While the full results of this search 
will be presented in a subsequent publication, here we 
focus on the dramatic outburst of HOPS 383. 

HOPS 383, at a = 5 h 35 m 29 s .81, S = -4°59'51"1 
(J2000), was identified as a protostar by Mcacat h et al.1 
42011 and is a source targeted by HOPS, the Her- 
schel Orion Proto star Survey (e.g., iStutz et all 120131 
iManoi et al .1120 1. 'll) . We ass ume it is at a d i stanc e of 420 
pc, the same as adopted bv lMeeeath et al.l (120121) for the 
entire Orion complex. 

2.1. Pre-Outburst Observations 

On 2000 December 12. iPeterson et, alJ (|2008H obtained 
a AT-band image of HOPS 383 with SQIID, the Simul¬ 
taneous Quad Infrared Imaging Device at the Kitt Peak 
2.1 m telescope. In 2004, HOPS 383 was detected at 4.5 
(im by the Infrared Array Camera (IRAC) and at 24 /jm 
by the Multiband Imaging Photometer (MIPS) aboard 
Spitzer. The IRAC observations were obtained in two 
epochs, one on 2004 March 9 and the other on 2004 Oc¬ 
tober 12, and the MIPS observation was on 2004 March 
20. An accounting of the Spitze r observations and sub ¬ 
seq uent analysis can be f ound in iKrvukova et all (|2012l) 
and iMeeeath et al.l (|2012f) . Color corrections for the pre- 
and post-outbust IRAC fluxes, which were determined 
for a flat-spectrum source, were estimated to be minimal. 
The pre- and post-outbust MIPS data were corrected by 
+5% relative to a Rayleigh-Jeans flux law due to the 
redness of the source. We also present a 450 pm image 
from SCUBA, the Submillimctre Commo n-User Bolome¬ 
ter Array, that was obtained in 1998 (i.Tohnstone fc Ballvl 

EM3). 

2.2. Post-Outburst Observations 

A /\ s -band image of HOPS 383 was obtained on 2009 
November 25 with NEWFIRM, the NOAO Extremely 
Wide-Field Infrared Imager, at the Kitt Peak 4 m tele¬ 
scope. A region including HOPS 383 was observed with 
the Spitzer InfraRed Spectrograph (IRS) in mapping 
mode on 2006 October 20. We performed photometry 
of HOPS 383 on the peak-up images acquired at 15.8 
and 22.3 pm; details will be provided in our forthcoming 
publication about the larger variability search. A sec¬ 
ond MIPS observation was acquired on 2008 April 19. 
As part of th e YSOVAR program (Young Stellar Ob ject 
VARiability; iMorales-Calderdn et alJl2011t lliebull et a. 1.1 
l20Ti . HOPS 383 was observed 81 times in the fall of 
2009 and 11 times in the fall of 2010 with the IRAC 3.6 
and 4.5 pm channels. 

WISE data were acquired in two visits on 2010 March 
8-9 and 2010 September 15-16. The first visit obtained 
photometry at 3.4, 4.6, and 22 (im and upper limits at 12 
pm. The second visit obtained photometry only at 3.4 
and 4.6 pm. We applied color corrections to the pho¬ 
tometry by fitting the model SED (§3.2) within each 
bandpass with a Planck function (bands 1 and 2) or a 
power law (bands 3 and 4) and interpolating between 
cases listed in the WISE Explanatory Supplement Q In¬ 
cluding the additional correction recommended at 22 pm, 
the corrections were —18%, —6.1%, +1.3%, and —12% in 
bands 1 though 4. 

12 See http://wise2.ipac.caltech.edu/docs/release/allsky/ 
expsup/sec4_4h. html. 


Table 1 

Photometry for HOPS 383 


A 

(pm) 

F u 

(mJy) 

cF u 

(mJy) 

Instrument 

Date 

2.2 

0.374 

0.0551 

SQIID 

2000 Dec 12 1 

2.2 

0.850 

0.0272 

NEWFIRM 

2009 Nov 25 1 

3.4 

0.585 

0.0969 

WISE 

2010 Mar 8-Sep 16 

3.6 

1.54 

0.0770 

IRAC 

2009 Oct 23 2 

4.5 

1.30 

0.0770 

IRAC 

2004 Mar 9-Oct 12 

4.5 

6.34 

0.317 

IRAC 

2009 Oct 23 2 

4.6 

10.5 

0.358 

WISE 

2010 Mar 8-Sep 16 

12 

< 16.5 


WISE 

2010 Mar 8-9 

16 

3.50 

0.450 

IRS 

2006 Oct 20 

22 

80.5 

3.32 

IRS 

2006 Oct 20 

22 

198 

16.8 

WISE 

2010 Mar 8-9 

24 

5.34 

0.282 

MIPS 

2004 Mar 20 

24 

187 

9.87 

MIPS 

2008 Apr 19 

70 

13100 

150 

PACS 

2010 Sep 10 2 

100 

25400 

1360 

PACS 

2010 Oct 8 

160 

35100 

2010 

PACS 

2010 Sep 10-28 

350 

6620 

2650 

SABOCA 

2011 Sep 16 

870 

996 

199 

LABOCA 

2010 Oct 24 


1 Flux in an aperture of radius 4.8" with a sky-subtraction annulus 
from 52" to 60". It likely contains a contribution from nebulosity 
not associated with HOPS 383; however, this should be constant 
between the two epochs. 

2 Representative of a time series. 

With Herschelwe observed HOPS 383 on 2010 Septem¬ 
ber 10 and 28 in the 70 pm and 160 pm bands of the 
Photodetector Array Camera and Spectrometer (PACS). 
Thes e observations are Groups 19 and 135 in Table 
2 of IStutz etahl ( 2013 0, who give additional informa¬ 
tion about the observations and describe the data pro¬ 
cessing and aperture photometry. We also include 100 
pm pho tometry from the Herschel Gould Belt Survey 
(jAndre et al .1120 1 ft observation of 2010 October 8 as well 
as 18 Herschel/PACS 70 pm data points obtai ned be¬ 
tween 2011 February 25 and 2012 August 27 (jBillot et al.1 

l20ia . 

For submillimeter coverage, we report beam fluxes 
from our Atacama Pathfinder Experiment (APEX) sur¬ 
vey of Orion. We acquired a 350 pm image on 2011 
September 16 with SABOCA, the Submillimetre APEX 
Bolometer Camera, and an 870 pm image on 2010 Octo¬ 
ber 24 with LABOCA, the Large APEX Bolometer Cam¬ 
era. Details of the APEX obser vations, dat a redu ction, 
and photometry are reported in IStutz et al.1 (l2013lh 

3. RESULTS 

3.1. Detection and Timing Constraints 

Fluxes for HOPS 383 before and after the outburst 
appear in Table [T] while Figure [1] shows near- and mid- 
infrared images. In the K band, the point source is not 
detected at either epoch, but nebulosity appears post¬ 
outburst. We report fluxes inside an aperture of radius 
4.8", the largest radius that does not include significant 
emission from nearby stars, with subtraction of the me¬ 
dian sky signal in an annulus extending from 52" to 60" 
that is free from bright point sources or extended emis¬ 
sion. These fluxes may include contamination from neb¬ 
ulosity in the region unrelated to HOPS 383, but this is 
likely to be constant across the two epochs, indicating an 
increase of 0.476 mJy at K due to the outburst. 

The best demonstration of the change in the point 
source appears in the MIPS images, where the source 
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Figure 1. Near- and mid-IR images of HOPS 383 before and after its luminosity increase. Top: Pre-outburst imaging from SQIID and 
Spitzer. Bottom: Post-outburst imaging from NEWFIRM, WISE , and Spitzer. The position of HOPS 383 is marked in green. 
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Figure 2. Submillimeter images. A SCUBA image at 450 pm 
obtained in 1998 appears on the left, a SABOCA image at 350 
pm obtained in 2011 appears in the center, and the ratio of the 
post-outburst to the pre-outburst image appears on the right. The 
position of HOPS 383 is marked in green. 

became a factor of 35 brighter at 24 /jm between 2004 
and 2008. The 2006 IRS measurement at 22 fin i is much 
closer to the 2008 MIPS flux density than the 2004 MIPS 
flux density; therefore, we conclude the outburst began 
between 2004 October 12 and 2006 October 20, with a 
subsequent rise in luminosity up to 2008. 

Figure [2] shows a SCUBA 450 fim image from 1998, 
a SABOCA 350 fim image from 2011, and the ratio of 
SABOCA to SCUBA. They show the appearance by 2011 
of a bright source at the position of HOPS 383. While 
the extended emission has a ratio of ~ 2.1 due to the 
difference in wavelength, the source has a ratio of ~ 4, 
implying an increase in brightness of > 2. 

The SED of HOPS 383 appears in Figure [3] Due to 
the paucity of pre-outburst data, we focus on the post¬ 
outburst SED. For data products where time series exist, 
the variability within the time series is small compared to 



Figure 3. Pre- and post-outburst photometry for HOPS 383, with 
pre-outburst data in blue and post-outburst data in red. The best- 
fit SED from the HOPS grid of radiative transfer models is shown 
with a gray curve. The WISE point with an arrow below it is an 
upper limit. For data products where time series exist, a represen¬ 
tative flux is shown. Green circles mark the data used in modeling. 

the range of the whole SED, so we choose a representative 
observation for the SED; the time series are discussed in 
Section 3.3. 

3.2. Properties of the Post-Outburst Source 

The bolometric lumin osity Lhni of the post-outburst 
SED is 7.5 Lq. iDunham et all (|2008fl showed that, 
for protostars, the correlation between flux density at a 
given wavelength and luminosity tightens with increasing 
wavelength between 3 and 70 fim. Lacking pre-outburst 
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70 pm data, we estimate that the change in luminosity 
was equal to the factor of 35 increase in the 24 /an flux 
density, giving a pre-outburst luminosity of 0.2 Lq. 

A common diagnostic of evolutionary state is the ra¬ 
tio of submillimeter luminosity L smm to bolometric lu¬ 
minosity. The Cla ss 0 objects have L sinm /Lb 0 \ > 0.5% 
(iAndre et al.lll993l) . and HOPS 383 has a ratio of 1.4%, 
confirming its deeply embedded status and making it 
the only known Class 0 accretion outburst. This ratio 
is toward the low end of the range observed for the ex- 
treme C lass 0 obje cts known as PACS Bright Red Sources 
( Stutz et al. 2013 :. 

Another diagnostic is the bolometric temperature, the 
temperature of a blackbody with the s ame mean fre¬ 
quency as that of the protostellar SED dMvers fe Laddl 
119931) . The Cla ss 0 protostars have Tboi < 70 K 
dChen et al.lll995l ). and these correspond roughly to pro- 
tostars in which the majority of the mass is in the 
infall ing envelope, not yet in the star (iDunham et al.l 
IM1 . The SED of HOPS 383 has 71,oi — 43 K, again 
consistent with Class 0. 

Additional evidence for a deeply embedded object 
comes from fitting a modi fied bla c kbody to the SED, 
where the method of iStutz et all ( 201 31 gives a peak 
wavelength of 106 /mi (larger peak wavelengths imply 
denser envelopes) and a lower limit to the envelope 
mass of 0.2 M 0 . Estimating the mass from the 870 
/mi APEX flux alone, assuming a temperature of 18 K 
(Stutz et al.. in preparation ) and OH5 opacities from 
lOssenkonf fo Hennind ( 19941 1. gives a mass of 0.7 M 0 in 
the LABOCA beam, which has a half-width at half max¬ 
imum corresponding to a radius of 4000 AU. These mass 
estimates point to a protostar in the earlier phases of 
infall. 

To estimate the luminosity, inclination angle, and cav¬ 
ity opening angle of HOPS 383, we modeled the Spitzer, 
Herschel , and APEX phot ometry of the source w ith the 
radiative transfer code of iWhitnev et al.1 (|2003H . The 
HOPS team created a grid of 3040 model SEDs with pa¬ 
rameters approp riate for protostars, first described by 
lAli et al.l (2010 1. and updated by E. Furlan et al. (in 
preparation). We found the best fit by minimizing R, 
which measures the logarithmic deviation of the models 
from t he observation s in units of the fractional uncer¬ 
tainty (iFischer et al.1120121) . 

The best-fit model (R = 1.7) has a total luminosity of 
8.7 Lq, an inclination of 41° from pole-on, and a cavity 
opening angle of 25°. Other models that provide sat¬ 
isfying fits to the data (i? < 2.1) have total luminosi¬ 
ties ranging from 6 to 14 Lq, inclinations ranging from 
41° to 63° from pole-on, and cavity opening angles from 
15° to 35°. (The total luminosity can differ from the 
observed luminosity due to the non-isotropic radiation 
field.) The scattered-light cone extending to the north¬ 
west and bright features to the southeast in the NEW- 
FIRM image of Figure [1] appear to be from radiation 
escaping the outflow cavity and imply an intermediate 
inclination angle, consistent with the SED. 

3.3. Post-Outburst Variability 

Our best sampling of the variability of HOPS 383 
comes from the Spitzer YSOVAR data at 3.6 and 4.5 /im 
and the Herschel/PACS 70 /im data. Figure [4] shows 
these light curves. After the 2 mag jump at 4.5 /im 
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Figure 4. Time-series photometry. Top two rows: Magnitude ver¬ 
sus time for the YSOVAR 2009 (left) and 2010 (right) campaigns. 
Bottom row : Flux density versus time for the Herschel /PACS 70 
pm data. The gray rectangle marks the epoch of the 2010 YSO¬ 
VAR campaign, and the two small panels zoom in on the two 2011 
PACS epochs. 


between 2004 and 2009 (not shown), the 2009 season 
of YSOVAR data yielded remarkably cons tant photom¬ 
etry, f ailing the variability tests laid out in IRebull et al.1 
(|2014j) . The slopes of the best-fit lines to the light curves 
indicate brightenings of only 0.07 mag at 3.6 pm. and 
0.08 mag at 4.5 pm, and the [3.6] — [4.5] color reddened 
by only 0.01 mag. 

There were two day-long series of WISE observations 
at 3.4 and 4.6 pm between the 2009 and 2010 YSO¬ 
VAR campaigns, on 2010 March 8-9 and 2010 September 
15-16. No significant variability from the WISE fluxes 
shown in the SED (Fig. [3]) was detected. The second 
(2010) YSOVAR season showed a slight fading of the 
source, with a dimming of 0.20 mag at 3.6 pm and 0.15 
mag at 4.5 pm. Again, the color was more consistent 
than the magnitudes, with a reddening of 0.05 mag over 
the window. 

The Herschel /PACS 70 pm data show variability at the 
18% level. The first significant gap in the PACS obser¬ 
vations coincides with the 2010 YSOVAR window; both 
light curves suggest a decline over this period. Between 
the first 70 pm observation on 2010 September 10 and 
the second-to-last one on 2012 March 18, the flux density 
declined by 13% before recovering on 2012 August 27 to 
its brightest yet. 


4. DISCUSSION 

We have presented multi-epoch infrared data from 2.2 
to 870 pm to show an outburst of the Class 0 protostar 
HOPS 383, adjacent to the NGC 1977 nebula in Orion, 
between 2004 and 2006. By 2008, the source became 35 
times brighter at 24 pm than it was in 2004. The post- 
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outburst luminosity, which is uncertain due to modeling 
degeneracies, is most likely between 6 and 14 Lq. Sub¬ 
sequent monitoring of the source at 3.6, 4.5, and 70 /im 
finds variability but no evidence for a significant decline 
in the source luminosity between 2009 and 2012. 

One possible mechanism for the observed rise in the 
24 (im flux density would be the removal from the line 
of sight of a large amount of extinguishing material. We 
are able to rule this o ut on two grounds. First, using 
th e opacity law from Ormell (120 1 II) plotted in Figure 3 
of lFischer et all ( 201 il l and a gas-to-dust ratio of 100, a 
factor of 35 increase in the 24 p ,m flux density would cor¬ 
respond to a reduction in the column density of 1.3 x 10 23 
cm -2 , which would correspond to an implausible reduc¬ 
tion in Ay of 70. Second, the source became brighter at 
submillimeter wavelengths (Fig. O, which cannot be the 
result of a drop in extinction and is most likely due to 
an incre ase in the envelope tem perature due to increased 
heating (iJohnstone et al.l [20131) . Ruling out a change in 
extinction implies that the observed flux increase was due 
to a genuine increase in the source luminosity. 

YSOs undergoing such luminosity outbursts are gen¬ 
erally classified as EX Lup objects, with repeated short¬ 
term flux increases of about a magnitude that persist for 
approximately one year, or as FU Ori objects, with a 
single burst of several magnitu des and the persistence o f 
an elevated state for decades dReipurth fe Asuinl 1 201 Of) . 
In both cases, an increase in the accretion rate onto the 
star is thought to be responsible for the luminosity in¬ 
crease. Outbursts have also been observed with light 
curves that do not fit cl eanly into either c ategory; e.g., 
V1647 Ori (lAspinl 12011) and V2492 Cyg (iCovev et all 
120111: iHillenbrand et ahf 20131! , and final confirmation of 
the type of outburst requires optical or near-infrared 
spectra to diagnose conditions in t he re gion w here disk 
material is accreting onto the star (iConnellev fc Oreend 
120101 !. Due to the lack of such spectra for a deeply em¬ 
bedded source, we are limited in our ability to conclu¬ 
sively identify the class of outburst or the precise mech¬ 
anism for the luminosity increase in HOPS 383, but the 
post-outburst light curves indicate the persistence of el¬ 
evated luminosity from the first evidence of brightening 
in 2006 to the most recent Herschel/ PACS imaging in 
2012, inconsistent with the short-term EX Lup events. 

HOPS 383 is unambiguously a Class 0 YSO based on 
the arguments presented in Section 3.2. While outburst- 
ing YSOs are known to have envelo pes in some cases 
(jQuanz et ;d.i 2007 : : Qreen et al.ll2013l ) , none has yet been 
detected that is as deeply embedded as HOPS 383. De¬ 
tections of d isks in such deep l y em bedded protostars 
remain rare dTobin et al.1 12012L 120131! . but there is in¬ 
direct evidence from the ubiquity of outflows and jets 
that there are di sks around nearly all protostars (e.g., 
iFrank et al.ll2014 j. Since proposed mechanisms for ac¬ 
cretion outbursts are gr o unded in disk in stabilities (e.g., 
iVorobvov fe Basul 120101 ; iZhu et all 120101! , the outburst 
of HOPS 383 is another form of indirect evidence for a 
disk in a Class 0 source. The discovery of this outburst 
demonstrates that episodic accretion can occur very early 
in the star formation process, and we encourage follow¬ 
up observations to begin to understand the physics at 
work. 
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